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ABSTRACT 

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole 
(NEP) deep field using recently-obtained, wider CEHT optical/near-IR images. AKARI 
has obtained deep images in the mid-infrared (IR), covering 0.6 deg^ of the NEP deep 
field. However, our previous work was limited to the central area of 0.25 deg^ due to 
the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, 
we recently obtained CEHT optical and near-IR images over the entire AKARI NEP 
deep field. These new CEHT images are used to derive accurate photometric redshifts, 
allowing us to fully exploit the whole AKARI NEP deep field. 

AKAREs deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 
7, 9, 11, 15, 18, and 24/im) exists nowhere else, due to filter gaps of other space tele¬ 
scopes. It allows us to estimate restframe 8jim and 12/im luminosities without using a 
large extrapolation based on spectral energy distribution (SED) fitting, which was the 
largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained 
more reliably due to the superior filter coverage. The resulting restframe 8/im, 12/rm, 
and TIR EEs at 0.15 < z < 2.2 are consistent with previous works, but with reduced 
uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. 
In terms of cosmic infrared luminosity density (Hir), we found that the Hir evolves as 
OC (1 + z)4-2±0.4^ 


1 INTRODUCTION 


Studies of the extragalactic background suggest that at least half 
the luminous energy generated by stars has be en reprocessed 
into the infrared (IR) emission by dust grains i Laeache et al.l 
1 19991 ; IPuget et ^ 1 19961 : iFranceschini et kj l2008h . indicating 
that dust-obscured star formation (SF) was more important at 
higher redshifts than today. 


iTakeuchi et alj LOOSah reported that the IR-to-UV lumi¬ 
nosity density ratio, pL{dust)/P l{uv) evolves from 3.75 (z=0) 
to 15.1 by z=l .0, af t er a car eful treatment of the sample selection 
effects. lOoto et afl LoiObh suggested that total infrared (TIR) 
luminosity accounts for ~70% of total star formation rate (SFR) 
density at 2=0.25, and 90% by 2=1.3. These results highlight 
the importance of probing cosmic star formation activity at high 
redshift in the infrared bands. 
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Figure 1. The 8fim luminosity distributions of samples used to compute 
restframe 8/xm LFs. Luminosity unit is given as the logarithmic with 
respect to solar luminosity (Lq). 

In our previous work l lGoto et alj|2010bh . we tried to ad¬ 
dress the i ssue using AKARFs m id-infrared data in the NEP 
deep field jMatsuhara et alj [20061) . However, the analysis was 
limited to 0.25 deg^, where we had deep optical imaging from 
Suharu/Suprime-Cam (Wada et al. in preparation), i.e., only 40% 
of the area was used for the analysis. To improve the situation, 
we have recently obtained much wider optical/near-IR imaging 
of the entire AKARI NEP deep field using CEHT MegaCam and 
WIRCam. 

AKARFs continuous mid-IR filters allow us to estimate 
the MIR (mid-infrared)-luminosity without using a large k- 
correction based on SED models. In this way, we are able to 
eliminate the largest uncertainty in previous works. By taking 
advantage of the unique filters and the larger field coverage, 
in this work we present more accurate estimates of restframe 
8, 12^m and total infrared (TIR) LFs using the AKARI NEP- 
Deep data. We define the total IR luminosity as the integral 
over the rest-frame wavelength range 8-1000/rm, which approx¬ 
imates the total luminosity emitted by interstellar dust, free from 
contamination by starlight. Unless otherwise stated, we assume 
a cosmology with {h, fim, Ha) = (0.7, 0.3, 0.7). 


2 DATA & ANALYSIS 

2.1 Multi-wavelength data in the AKARI NEP Deep field 

To obtain t he TIR LFs, we adopt the method presented by 
iGoto et alJ ( l20I0lj) . Major improvements over our previous 
work can be summarized as follows: 

• The new AKARI NEP catalog. We have t he advantage 

of the new AKARI 9-band mid-IR catalog in iMurata et al.l 
(l2013h . which improved the corrections for scattered and stray 
light, artificial patterns, and artificial sources due to bright 
sources. As a result, we achieved deeper 5a detection limits, 
which are 11, 9, 10, 30, 34, 57, 87, 93, and 256 /rJy in the 
N2, N3, NA, S7, S9W, Sll, L15, L18W, and L24 filters, re¬ 
spectively. _ 

• Wider field coverage. Our previous work jGoto et al.l 


Table 1. Summary of filters used for LFs. 


Estimate 

Redshift 

Filter 

8/im LF 

0.28<z<0.47 

SlI (Il/rm) 

8/im LF 

0.65<z<0.90 

L15 {I5gm) 

Hfim LF 

1.09<z<1.41 

LI8W(18Aim) 

8/im LF 

L78<z<2.22 

L24 (24/im) 

12//m LF 

0.15<z<0.35 

LI5 (IS/rm) 

llfim LF 

0.38<z<0.62 

L18W(I8/im) 

\2fim LF 

0.84<z<1.16 

L24 (24/rm) 

TIRLF 

0.2<z<0.5 

S7, S9W, Sll, L15, L18VF, and L24 

TIRLF 

0.5<z<0.8 

S7, S9W, Sll, L15, L18VF, and L24 

TIRLF 

0.8<z<1.2 

S7, S9VF, Sll, L15, L18VF, and L24 

TIRLF 

1.2<z<1.6 

S7, S9W, Sll, L15, L18VF, and L24 


l2010bl) was limited to the central area of 0.25 deg^, where we 
had deep optical data, making the results vulnerable to cosmic 
variance. We recently obtained CFHT/MegaCam g, r, i, z, Y, J 
and, K photometry (PLGoto), which covers the entire AKARI 
NEP deep field (0.9 deg^). The wider coverage allows us to fully 
exploit the unique space-home mid-IR data, while reducing the 
uncertainty from cosmic variance. 

• More accurate photometric redshifts. Due to the availability 
of d eeper near-IR photo metry, combined with existing u-band 
data l lTakagi et al.ll201 j) . our photometric redshift accuracy was 
also impro ved to ^ 0.032 at z < 1, and ~0.117 at 
z > 1 (see lOi et alT2014l for details of photo-z calculation). 

Thanks to these improvements, we have a larger sample of 
5761 infrared galaxies. AKARFs 9 mid-IR filters are still ad¬ 
vantageous in providing us with a unique continuous coverage 
at 2-24^m, where there is a gap between the Spitzer IRAC and 
MIPS, and the ISO LW2 and LW3. 


2.2 Analysis 

We compute LFs using the 1/Umax method, as we did in 
iGoto et ^ ( l2010bh . Uncertainties of the LF values stem from 
various factors such as fluctuations in the number of sources in 
each luminosity bin, the photometric redshift uncertainties, the 
fc-correction uncertainties, and the flux errors. To compute the 
errors of LFs we performed Monte Carlo simulations by cre¬ 
ating 1000 simulated catalogs. Each simulated catalog contains 
the same number of sources, but we assigned a new redshift to 
each source, by following a Gaussian distribution centered at the 
photometric redshift with the measured dispersion A«/(l + z). 
The flux of each source is also changed; the new fluxes vary 
according to the measured flux error following a Gaussian dis¬ 
tribution. 

For the 8^m and the 12/im LFs, we can ignore the errors 
due to the fc-correction thanks to the continuous AKARI MIR 
filter coverage. The TIR errors are estimated by re-performing 
the SED fitting for each of the 1000 simulated catalogs (see be¬ 
low for more details of the SED fitting). 

We did not consider the uncertainty related to cosmic vari¬ 
ance here since our field coverage has been significantly im¬ 
proved. For our analy sis, each redshif t bin c overs 4 x 10® 
Mpc® of volume. See lMatsuhara et ^ ( l2006l) for more discus¬ 
sion on the cosmic variance in the NEP field. 

All the other errors described above are added to the Pois¬ 
son errors for each LF bin in quadrature. 
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3 RESULTS 


and Ltir as 


3.1 The 8^im LF 

Monochromatic 8^m luminosity (Lsu m) is known to cor¬ 
relate well v yith t h e TIR lurn i nosity teabhedge et alj 1 20061 : 
iHuang et"^ IiOOTI : lOoto et aP 1201 lah . especially for star¬ 
forming galaxies, because the rest-frame 8^m flux is dominated 
by prominent PAH (polycyclic aromatic hydrocarbon) features 
such as those at 6.2, 7.7, and 8.6 ^m. 

Since AKARI has continuous coverage in the mid-IR wave¬ 
length range, the restframe 8/rm luminosity can be obtained 
without a large uncertainty in fc-correction at the corresponding 
redshift and filter. For example, at 2=0.375, restframe 8^m is 
redshifted into Sll filter. Similarly, LIS, LISFF and L24 cover 
restframe 8^m at 2=0.775, 1.25 and 2. This filter coverage is an 
advantage with AKARI data. Often in previous work, SED mod¬ 
els were used to extrapolate from Spitzer 24/rm fluxes, produc¬ 
ing the largest uncertainty. This is not the case for the analysis 
present in this paper. Table [T] summarizes all filters used for our 
computation. 

To obtain the restframe 8/im LF, we used sources down 
to 80% completen ess limits in each band as measured in 
iMurata et alj J 2 OI 3 I) . We excluded those galaxies whose SEDs 
are better fit with QSO templates. This removed 2% of galaxies 
from the sample. _ 

W e used the completeness curve presented in iMurata et al.l 
( l2013h to correct for the incompleteness of the detections. How¬ 
ever, this correction is 25% at maximum, since our sample is 
brighter than the 80% completeness limits. Our main conclu¬ 
sions are not affected by this incompleteness correction. To 
compensate for the increasing uncertainty at increasing 2 , we 
use four redshift bins of 0.28< 2 <0.47, 0.65< 2 <0.90, 
1.09< 2 <1.41, and 1.78< 2 <2.22. We show the Lsum distri¬ 
bution in each redshift bin in Fig.[T] Within each redshift bin, we 
use the 1/Umax method to compensate for the flux limit in each 
filter. 

We show the computed restframe 8/rm LF in Fig. The 
an'ows mark the 8/im luminosity corresponding to the flux limit 
at the central redshift in each redshift bin. Error bars on each 
point are based on the Monte Carl o simulation (S 12.2b . and are 
smaller than in our previous work dOoto et al.ll2010blL To com¬ 
pare with previous work, the dark-yellow dot-dashed line also 
sh ows the 8/rm LF of star-forming galaxies at 0 < 2 < 0.3 
bv iHuang et al.l d2007l) . using the 1/Unax method applied to the 
IRAC 8^m GTO data. Compared to the local LF, our 8^m LFs 
show strong evolution in luminosity . At higher redsh i ft, our 8/im 
LFs are located be tween those by ISabbedge et al.l ( l2006h and 
ICaputi et alj d2007h. At z=0.7, our LF is in a perfect agreement 


Ltir = (20 ± 5) X t/L°'g^*° °^(±44%). 


( 1 ) 


with iFu et al.l d2010f) . which is based on Spitzer IRS spectra. 


3.2 Total IR luminosity density based on the S/rm LF 

One of the important purposes in computing IR LFs is to esti¬ 
mate the IR luminosity density, which in turn is an extinction- 
free e stimator of the cosmic star formation density dKennicuttl 
Il998h . 

We estimate the total infrared luminosity density by inte¬ 
grating the LF weighted by the luminosity. First, we need to 
convert Ls^m to the total infrared luminosity. 

A strong correlation between Ls^m and tota l infrared lu¬ 
minosity (Z/tir) has been r eported in the literature dCaputi et al.l 
l2007l : l^vouzet et ^l2008h . Using a large sample of 605 galax- 
ies detected in the far-infrared by the AKARI all sky survey, 
lGotoetalJd2011bl) estimated the best-fit relation between Lsfj,m 


The Ltir is based on AKARFs far-IR photometry in 
65,90,140, and 160 /rm, and the measurement is based 

on AKARFs 9fim photometry. Given the superior statistics and 
availability over longer wavelengths (140 and 160/rm), we used 
this equation to convert into Ltir- 

However, the above-mentioned conversion is based on 
local star-forming galaxie s. It is uncertain wh ether it holds 
at higher redshift or not. iNordon et al.l d2012h reported that 
“main-sequence” galaxies (normal star-forming galaxies that 
follow the stellar-mass and SFR relation) tend to have a sim¬ 
ilar I/g^m/LTiR regardless of Lir and redshift, up to z~2.5, 
and Lgfim/LTiR decreases with increasing offset above the main 
sequence. They suggested that this is due to a change in the 
ratio of polycy clic aromatic hydroca rbon (PAH) to Ltir- On 
the other hand, ISavouzet et alj ( l2008h checked this by stacking 
24/im sources at 1.3 < 2 < 2.3 in the GOODS fields to find 
whether the stacked sources are consistent with the local rela¬ 
tion. They concluded that the correlation is valid to link Lsfj,m. 
and I/TiR at 1.3 < 2 < 2.3. lTakagiet'ZI 1 I 2 OIOI) also showed 
that local Lr,rfj,m vs Lti r relation holds t rue for IR galaxies 
at z~l (see their Fig. 10). fPope et af] ( l2008h showed that 2 ~2 
sub-millimeter galaxies lie on the relation between Ltir and 
LpAH , 7.7 that has been established for l ocal starburst galaxies . 
The S'7o/<S'24 ratios of 70/rm sources in lPapovich et alj ( l2007h 
are also consistent with the local SED templates. Following 
these results, we use equation (D for our sample. However, we 
should keep in mind th at there may be a p o ssible evolution wit h 
redshift as discussed in iRigbv et al.l l l200^ : lHuang et all ( l2009h . 

The conversion, however, has been the largest source of 
error in estimating Ltir v alues from Lsum - Reported disper¬ 
sions a re 37 , 55 and 44% bviBavouzet et all ( l2008h . lCaputi et all 
l l2007l) . and iGoto et all l l2011bh . respectively. It should be kept 
in mind that the restframe 8/im is sensitive to the star-formation 
activity, but at the same time, it is where the SED models have 
strongest discrepancies due to the complicated PAH emission 
lines. Possible SED evolution, and the presence of (unremoved) 
AGN will induce further uncertainty. A detailed comparison 
of diff erent conversions is presented in Fig.12 of ICaputi et alj 
l l2007l) . who reported a factor of -^5 differences among various 
models. 

Having these cautions in mind, the 8/rm LF is weighted 
by the Ltir and integrated to obtain the TIR density. For in¬ 
tegration, we first fitted an analytical function to the LFs. In 
the literature, IR LFs were fit better by a do uble-power law 
jBabbed^e_et^ 2006h or a double-exponentia l jSaunders_et^ 
19901: IPozzi et al.l 12004 iTakeuchi et all l2006l: iLe Floc’h et ~ 


2005h than a Schechter function, which declines too steeply at 
high luminosities, underestimating the number of bright galax- 
ies. In this work, we fi t the 8/im LFs using a double-power law 
iBabbedge et al.|[2006l) as is shown below. 


^(L)dL/L* = ^ 

' L* 


$(L)dL/L* = 


dL/L*, {L < L*) 


1-/3 


dL/L*, {L > L*) 


( 2 ) 


(3) 


First, the double-power law is fitted to the lowest redshift LF at 
0.28< 2 <0.47 to determine the normalization (4>*) and slopes 
(a, P). For higher redshifts we do not have enough statistics 
to simultaneously fit 4 parameters ($*, L*, a, and /3). There¬ 
fore, we fixed the slopes and normalization at the local values 
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0.28<z<0,47 


0.65<z<0.90 




1 .09<z<l .41 



1.78<z<2.22 



Figure 2. Restframe 8/i.m LFs based on the AKARI NEP-Deep field. The top panel shows results from all the redshift bins. The succeeding panels 
show them one by one for clarity. The blue diamonds, the purple triangles, the red squares, and the orange crosses show the 8/i.m LFs at 0.28 < z < 
0.47, 0.65 < 2 : < 0.90,1.09 < 2 : < 1.41, and 1.78 < 2 : < 2.22, respectively. AKARFs MIR filters can observe restframe 8//m at these redshifts 
in a corresponding filter. Error bars are estimated from the Monte Caro simulations ( ^2.21 . The dotted lines show analytical fits with a double-power 
law. V ertical arrows show the 8/im lu minosity correspond ing to the flux limit at the centr al redshift in each re dshift bin. Overplotted are lSabbedge et all 
(lOPj) i n the p ink dash-dotted lines. ICanuti et alJ i20Q7h in the cyan dash-dotted lines. iHuang et al.l j2007h in the dark-yellow dash-dotted lines, and 
IFu et all f2ol4 in the green dash-dotted line. 


and varied only L* for the higher-redshift LFs. Fixing the faint- 
end slope is a co mmon procedure with the depth of curre nt IR 
satellite surveys dSabhedge et al.ll2006l : ICaDuti et alj|2007h . The 
stronger evo lution in luminosity than in density found by pre- 
vious work jPerez-Gonzalez et al.ll2005l : iLe Floc’h et aljFoOSll 
also justifies this parametrization. The best-fit parameters are 
presented in Table We found that our evolves as oc 

(1 + in the range of 0.48 < z < 2. Once the best-fit 


parameters are found, we integrate the double power law outside 
the luminosity range in which we have data to estimate the TIR 
luminosity density, Otir- 


The hnal resulting total luminosity density (Oir) is shown 
in i l3.6l along with results obtained from other wavelengths. 
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Table 2. Best fit parameters for 8,12/im and TIR LFs 


Redshift 

LF 

L* ( 101 “Lq) 

<t>*(10“®Mpc“®dex“l) 

a 

/3 

0.28<z<0.47 

8/im 

ry O + 0.2 
^•'^-0.2 

^•^- 0.3 

1 CT+O.Ol 
i-'^-J_0.02 

0 4+0.6 

'^■^- 0.5 

0.65<z<0.90 

8/im 

0 1+0.1 

1.6 

1.53 

3.4 

L09<z<L41 

8/im 

S-Sloi 

1.6 

1.53 

3.4 

L78<z<2.22 

8/im 

5 . 7 ^ 5 :? 

1.6 

1.53 

3.4 

0.15<z<0.35 

12/im 

1 ]^+0.2 

1 1+0-2 
- 0.3 

1 gi+0.05 
^■^^-0.04 

3 4+1.4 

•^•^- 0.4 

0.38<z<0.62 

12/im 

7 4+0-2 
"^•^-0.2 

1.1 

1.91 

3.4 

0.84<z<L16 

12/im 

4 

^•°- 0.3 

1.1 

1.91 

3.4 

0.2<z<0.5 

Total 

loli 

T 4+0.2 

1 4 + 0 + 

0 O+0.2 
■^■■^-0.1 

0.5<z<0.8 

Total 

1 + 2 

3.4 

1.4 

3.3 

0.8<z<L2 

Total 

3+2 

3.4 

1.4 

3.3 

L2<z<L6 

Total 

1401“ 

3.4 

1.4 

3.3 



og [i/L^(12/xm) (Lq)] 


Figure 3. The ll/im luminosity distributions of samples used to com¬ 
pute restframe 12^m LFs. The Luminosity unit is in logarithmic solar 
units (Lq). 


3.3 12/im LF 

The 12/im luminosity (Li 2 ^m) has been well studied through 
ISO and I RAS. It is known to correlate closely with the T IR 
luminosity dSninoglio et alJl995l : IPerez-Gonzalez et al.ll2005h . 

As was the case for the 8/im LF, it is advantageous that 
AKARFs continuous filters in the mid-IR allow us to estimate 
restframe 12/im luminosity without much extrapolation based 
on SED models. 

At targeted redshifts of 2=0.25, 0.5, and 1, the LIS, L18VF 
and L24 filters cover the restframe 12/im, respectively. We sum¬ 
marize the filters used in Table [T] The methodology is the same 
as for the 8/im LF; we used the sample down to the 80% com¬ 
pleteness limit, corrected for the incompleteness, then used the 
1/Fmax method to compute the LF in each redshift bin. The 
Li 2 ^im distribution is presented in Fig.[3 The resulting 12/im LF 
is shown in Fig. [4] The light green dash-dotted line shows 12/im 
LF based on 89 3 galaxies at 0 < 2 < 0.3 in the IRAS Faint 
Source Catalog dRush et al.lll993h . The dark green dash-dotted 
line shows 12/im LF at 0.006 < z < 0.05 based on 223,982 
galaxies from WISE sources in Table 7 of iToba et all (l2014ll . 


Compared with these 2=0 LFs, the 12/im LFs show steady evo¬ 
lution with increasing redshift. In the range of 0.25 < 2 < 1, 
the Li 2 ;im evolves as oc (1 + 2 )^ ®^° '^. 

3.4 Total IR luminosity density based on the 12/im LF 

The 12/im is one of the most frequently used monochro¬ 
matic fluxes to estimate Ltir. The total infrared luminosity 
can be compu t ed from the Li 2 um using the co nversion in 
ICharv & Elbazl(l200lh : lpirez-Gonzalez et al.l ( l2005h . 

log Ltir = log(0.89lSi?) + 1.094 log Lia^^ (4) 

iTakeuchi et alj ( l2005bl) independently estimated the rela¬ 
tion to be 

log Ltir = 1.02 + 0.972 log Li2;im, (5) 

which we also use to check our Ltir — Li 2 ij,m conversion. 
As both authors state, these conversions contain an error of a 
factor of 2-3. Therefore, we should avoid conclusions that could 
be affected by such errors. The computed total luminosity den¬ 
sity based on the 12/im LF using equation (ID will be presented 
in i l3.6l 

3.5 TIR LF 

AKARFs continuous mid-IR coverage is also superior for SED- 
fitting to estimate Ltir. This is because for star-forming galax¬ 
ies, the mid-IR part of the IR SED is dominated by the PAH 
emiss ion lines, which reflect the SFR of galaxies dCenzel et alj 
Il998h . and thus, correlates well with Ltir, which is also a good 
indicator of the galaxy SFR. 

After photometric redshifts are estimated using the UV- 
optical-NIR ph otometry, we fix the redshift at the photo -2 
dOi et al.ll2014h . then use the same LePhare code to fit the 
infrared part of the SED to estimate TIR luminosity. We used 
iLagache et alj d2003lF s SED templates to fit the photometry us¬ 
ing the AKARI bands at >6/tm (57, S'9W, S'!!, L15, L18FF 
and L24). 

In the mid-IR, color-correction could be large when strong 
PAH emissions shift into the bandpass (a factor of ~3). How¬ 
ever, during the SED fitting, we integrate the flux over the band¬ 
pass weighted by the response function. Therefore, we do not 
use the flux at a fixed wavelength. As such, the color-correction 
is negligible in our process (a few percent at most). We also 
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0.15<z<0.35 




0.38<z<0,62 



0.84<z<l .1 6 



Figure 4. Restframe 12^m LFs based on the AKARI NEP-Deep field. Luminosity unit is logarithmic solar luminosity (Lq). The blue diamonds, the 
purple triangles, and the red squares show the 12/rm LFs at 0.15 < z < 0.35, 0.38 < z < 0.62, and 0.84 < z < 1.16, respectively. In the upper-left 
panel, we plot results from redshifts together. In the succeeding panels, we compare results from different redshifts separat ely for clarity. Vertical arrows 
show the 12/rm luminosity corresponding to the flux li mit at the central r edshift in each redshift bin. Overplotted are IPerez-Gonzalez et alj j2005h 
at 2=0.3, 0 . 5 and 0.9 in the dark-cyan dash-dotted lines. iToba et alj i2014li at 0< z <0.0 5 based on WISE in the dark green dash-do tted lines, and 
iRush et ^ il993h at 0< z <0.3 in the light green dash-dotted lines. Note lRush et alj jl993h is at higher redshifts than lXoba et all i2014h . 


checked that using d ifferent SED models dCharv & ElbazlIlOoi] : 
IPale & Heloijl2003l does not change our essential results. 


Table 4. TIR luminosity density values based on 8/zm LFs, shown with 
the brown stars in Fig.|6] 


Galaxies in the targeted redshift range are best sampled 
in t he 18p,m band due to the wide bandpass of the LISW fil¬ 
ter dMatsuhara et al.ll2006h . In fact, in a single-band detection, 
the 18^m image returns the largest number of sources. There¬ 
fore, we applied the 1/Vmax method using the detection limit at 
L18VK. We also checked that usin g the L15 flux limit d oes not 
change our main results. The same Lagache et al.l d2003h ’s mod¬ 
els are also used for fc-corrections necessary to compute I4nax 
and Knin-The redshift bins used are 0.2< z: <0.5,0.5< z <0.8, 
0.8< z <1.2, and 1.2< z <1.6. 


z ^tiRs^„{LqMpc 3/io8) 


0.48 

3.5± 1.7 

0.77 

4.8± 2.0 

1.25 

8.5±3.5 

2.00 

8.8±5.1 


Table 5. TIR luminosity density values based on 12/rm LFs, shown with 
the pink filled-triangles in Fig.|6] 


The obtained Ltir LFs are shown in Fig. [5] The uncertain¬ 
ties are estimated through the Monte Carlo simu l ations ( ^2.21 . 
For a local benchmark, we overplot iGoto et aP d2011ah . who 
cross-correlated the SDSS and the AKARI FIR all sky survey 
to derive local IR LFs with 2357 galaxies. Readers are also ref¬ 
ereed to Marchetti et al. (submitted) for local LFs based on Her- 
schel HerMES wide fields. The TIR LFs show a strong evolution 
compared to local LFs. At 0.25 < z < 1.3, Ltir evolves as 
oc(l + z)^'»±°-i. 

Using the same methodology as in previous sections, we 
integrate the Ltir LFs in Fig.j^through a double-power law fit 
(eq. |^and[3' The resulting total luminosity density (Gir) will 
be shown in 113.61 along with results from other wavelengths. 




0.25 

2.4± 1.1 

0.50 

5.7± 2.9 

1.00 

1L8±5.9 


3.6 Evolution of Gir 

In this section, we measure the evolution of Gir as a function 
of redshift. In Fig.[^ we plot Gir estimated from the TIR LFs 
(red circles, Table[3, 8/rm LFs (brown stars, Table|4j, and llptm 
LFs (pink filled triangles, Table|5j, as a function of redshift. Es¬ 
timates based on 12/im LFs and TIR LFs agree with each other 
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0.2<z<0,5 



0.5<z<0.8 



X 

CD 

TD 


0.8<z<1.2 


1.2<z<1.6 




Figure 5. The TIR LFs from the SED fit. The top panel shows results from all the redshift bins. The succeeding panels show one redshift bin at a time 
for clarity. Vertical arrows show the luminosity correspondin g to the flux limit at the central redshift in each redshift bin. We overplot z=Q IR LF ba sed 
on the AKARI FIR all sky surv ey in the black dash-do t line IGoto et ^|2Q1 lah . Overplotte d previous studies ar e taken from iLe Floc’h et"^ ^2QQ5t) in 
the dai'k-green, dash-dotted line. fMagnelli et ^ j2009h in the dai'k-red, dash-dotted line, and iHuvnh et alj i2QQ7h in the dark-yellow, dash-dotted line at 
several redshifts as mai'ked in the figure. 


Table 3. TIR luminosity density values and contributions by LIRG and ULIRG as a function redshift as in Fig.[^ 


z 

^Tm{LQMpc ®/10®) 

^lirg(LqMpc 3/10®) 

CiuLIRG{L@Mpc 3/10®) 

0.35 

3.7±2.1 

1.25±0.69 

0.07± 0.06 

0.65 

6.4± 3.0 

3.24±1.65 

0.23±0.13 

1.00 

10.9± 4.9 

6.61±3.01 

0.80± 0.39 

1.40 

22.0±11.3 

13.51±6.82 

3.97± 2.53 
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z 


Figure 6. The evolution of the TIR luminosity density based on TIR LFs (red circles, Table[3), 8/im LFs (stars. Tabled, and \2 ^ m LFs (filled tr i angles , 
Table[3. Values ai'e presented in Tabled Overplo t ted dot-dashed lines are estimates from the literat ure: iLe Floc’h et alj i20Q5l).lMagnelli et^ i2009l) . 
IPerez^onzalez et alJ ^20051) . ICaouti etalJ j2007h . iGrupoioni et alJ i2Q13h . iRodishiero et alJ j201Qt) . and iBabb ed^^ ^lJ j2QQ6h are in cy an, yellow, 
green, navy, dark green, orange, and pink, respectively. The purple dash-dotted line sh ows the UV estimate bv ISchiminovich et alj i2Q05h. The pin k 
dashed line shows t he total estimate o f IR (TIR LF) and UV Jschiminovich et ^l2QQ5l) . The open triangles are low-z results from lGoto et Si bOl lah : 
iHuang et alJi20Q7h : lToba et alJi2Ql4) in TIR, 8//m and 12//m, respectively. 



Z 


Figure 7. The evolution of the TIR luminosity density based on TIR LFs (red circles. Table [3). The blue open squares and orange filled squares 
corresponds to LIRG and ULIRGs samples only, also based on our -Lxir LFs. Plotted values are presented in Table[3] 
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very well, while those from 8 ^m LFs are slightly lower by a 
factor of a few than others. However, they are within one sigma 
of each other. Our 8 Atm LFs flatten s at 2 > 1.4. This is in good 
agreement with lCanuti et al.l l l2007h as shown with the dark blue, 
dash-dotted line. 

We compare our resu l ts with various papers in the litera- 
ture ^Le^locTi_et_aL[j 200 ^jBabbedeeet^i 2 l 200 ^JCanutie^^ 


2007|^ Perez;Gonzalezet^ _ 2005; _ Magnelli et al. 20091 : 

Gruppioni et al.l 20131 : Rodighiero et 7 iJ I2010I) in the dash- 
dotted lines. Our Oih, shows a very good agreement with these 
at 0 < 2 < 1 . 2 , as almost all of the dash-dotted lines lying 
within our error bars of Oir, from Ltir and 12/rm LFs. This is 
because an estimate of an integrated value such as Oir is more 
reliable than that of LFs. 

We parametrize the evolution of Oir using the following 
function. 


Oir(2) OC (1 -I- zf 


(6) 


By fitting this to the Oir, we obtained 7 = 4.2zfc 0.4. This is con¬ 
sistent with most previous works. For example, LeHoc^het^ 
( l2005h obtained 7 = 3.9 ± 0.4, IPerez-Gonzalez et al.l l l2005l) 
obtained 7 = 4.0 zb 0.2, Isabbedge et al. l l2006tl obtained 


7 = 4.5jln'7 ^fa giisUi et alJ l l2009ll obtained 7 = 3.6 ± 0.4. 
iGruppioni et alTlT OlSh: iBurg^ella et al.l ( 1201 3h obtained 7 = 


3.55 ± 0.1. Rodighiero et al.l l 201(lh obtained 7 = 3.55 ± 0.1. 
The agreement confirms a strong evolution of Gir. 


3.7 Differential evolution among ULIRG and LIRG 

In Fig. [7] we plot the contributions to Hir from LIRGs (Lu¬ 
minous InfraRed Galaxies; Ltir > IO^^Lq) and ULIRGs 
(Ultra-Luminous InfraRed Galaxies: Ltir > IO^^Lq, mea¬ 
sured from TIR LFs), with the blue open squares and orange 
filled squares, respectively. Both LIRGs and ULIRGs show 
strong evolution. From 2 = 0.35 to 2 = 1 . 4 , Hir by LIRGs in¬ 
creases by a factor of ~L 8 , and Hir by ULIRGs increases 
by a factor of ~10. The physical origin of ULIRGs in the lo¬ 
cal Universe is often a merger/ i nteraction dSanders & Mirabell 
Il996l : l^niguchi & Shiova|[T998l : lGotoll2005h . It would be inter¬ 

esting to investigate whether the merger rate also increases in 
proportion to the ULIRG fraction, or if different mechanisms 


can also produce ULIRGs at higher redshift (see early work by 

Swinbank et al. 

201(]l:lKartalteDe et al.l2012l:lTargett et al.l2013l: 

Stott etal.ll2013 

). 


3.8 Comparison to the UV estimate 

We have been emphasizing the importance of IR probes of the 
total SFR density (SFRD) of the Universe. However, the IR es¬ 
timates do not take into account the contribution of the unab¬ 
sorbed UV light produced by young stars. Therefore, it is im¬ 
portant to estimate how significant this direct UV contribution 

is. _ 

ISchiminovich et ai] l l2005tl based on the GALEX data sup¬ 
plemented by the VVDS spectroscopic redshifts, found that the 
energy density measured at 1500A evolves as oc (1 -I- 
at 0 < 2 < 1 and oc (1 -I- at 2 > 1. We overplot 

their UV estimate of psFR with the purple dot-dashed line in 
Fig.| 6 ] The UV estimate is almost a factor of 10 smaller than 
the IR estimate at most of the redshifts, confirming the impor¬ 
tance of IR probes when investigating the evolution of the total 
cosmic star formation density. In Fig. we also plot the total 
SFD (or Qtotai ) by adding Guv and Otir, with the magenta 
dashed line. In Fig.[ 8 l we show the ratio of the IR-to-UV SFRD 



Figure 8 . Htir to Duv ratio shown as a function of redshift. 


(Htir/ Guv) as a func tion of redshift. Althou gh the errors are 
large, Fig.[ 8 ]agrees with T^euchi_et_^ (1200^ and recent work 


based on Herschel bv iBurgarella et al 


ll2013h . Our result 


gests that Htir/ Guv evolves strongly with redshift: from 5.3 
at z=0.35 to 13.5 at z=L4, i.e., Htir explains already 85% of 
t^totai{= Huv + Utir) at 2=0.25, and that by 2 = 1 . 4 , 93% of 
the cosmic SFD is explained by the infrared. This implies that 
Htir provides a good approximation of the iltotai at 2 > 1 . 


4 DISCUSSION 

In this work, we have expanded our previous work on mid- 
IR LFs to the entire AKARI NEP deep field, using our new 
wide coverage CEHT optical/near-IR photometry. Our LFs are 
in good agre ement with our pre vious work based on a smaller 
(40%) area dOoto et al.l 1201 Obi) , but with much reduced er¬ 
rors. Our previous work had notable shot noise, which made 
it difficult to disc uss evolutionary and environmental effects 
dOoto et alj|2010ah . The LFs in this work are more stable, and 
allow us more precise discussion of the redshift evolution. This 
highlights the importance of wide field coverage, not only to 
reduce statistical errors, but also to overcome cosmic variance. 
AKARI has also obtained sha llower but wider NEP wide-field 
data, which covered 5.4 deg^ dLee et al ]|2009l: I Jeon et alj|2014l: 
iKim et aDl2012l) . We are making an effort to probe the bright- 
end of the LFs from these data, using spectroscopic redshifts 
dKimll201^ . In parallel, we are trying to obtain accurate photo¬ 
metric redshifts over the NEP-wide field. Initial data in r-band 
have just been taken using the new Subaru Hyper-Suprime Cam 
(r ~27mag, PI: Goto). 

By combining Herschel/SPIRE far-IR data, more accurate 
estimate of Ltir can be obtained than the mid-IR only esti¬ 
mation (e.g., iKarouzos et alJl2oT3) . In the AKARI NEP deep 
field, there exists Herschel data (PI: Sergeant). However, SPIRE 
data in the field is relatively shallow. The SPIRE source density 
is only 10% of that of AKARI sources. Only 1.1% of AKARI 
sources had a counterpart in SPIRE in a simple positional match 
within 5”. Therefore, if we want to exploit AKARFs data to the 
faintest flux limit, we need to proceed without SPIRE detection 
for 99% of the sources. Alternately, if we limit our sample to 
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SPIRE detected sources, we can only use 1.1% of the AKARI 
sources. One of the main purposes of obtaining the new CFHT 
imaging data was to increase the statistics. Therefore, in this pa¬ 
per, we decided not to limit our sample to Herschel detected 
sources. 

In other de eper fields, however, Herschel has revealed th e 
IR LFs to z=4 dGruppioni et 'ZI I 2 OI 3 I : lBurgarellaetZll2013h . 
This achievement has changed the role of the mid-IR studies; 
before Herschel, the mid-IR analysis was an effective way to es¬ 
timate SFR in the high-z Universe, even considering the large 
uncertainty in converting mid-IR photometry to Ltir or SFR. 
However, now that Herschel has provided direct far-IR measure¬ 
ments, the importance of mid-IR data has shifted toward under¬ 
standing more details of PAH or hot-dust emission, in compar¬ 
ison with far-IR data, in order to understand the SED evolu¬ 
tion be tween the mid-IR and far-IR. Many such attempts have 
begun l Elbazet ^ l201ll; iNordon et al.ll2012l ; ISuat et alj[20l^; 
iMurata et alj|2014 ). This work provides an important step from 
the mid-IR side, based on AKARFs unique continuous 9-band 
mid-IR coverage over a wide area (0.6 deg^). 


5 SUMMARY 

We have used recently obtained wide-held optical/near-IR imag¬ 
ing from CFHT. These data covered the entire AKARI NEP deep 
held, allowing us to utilize, for the hrst time, all of the infrared 
data AKARI has taken in the held. 

AKARFs advantage is its continuous hlter coverage in the 
mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24pm). 
No other telescope including Spitzer and WISE has such a hlter 
coverage in the mid-infrared. The data allowed us to estimate 
mid-IR luminosity without a large extrapolation based on SED 
models, which was the largest uncertainty in previous studies. 
Even for Ltir, the SED htting is more reliable due to this con¬ 
tinuous coverage of mid-IR wavelengths. 

We have estimated restframe 8pm, 12pm, and total infrared 
luminosity functions, as functions of redshift. All LFs are in 
agreement once converted to a TIR LF, and show strong evolu¬ 
tion from 2=0.35 to 2=2.2. Dir evolves as oc {1 + We 

found that the ULIRG (LIRG) contribution increases by a fac¬ 
tor of 10 (1.8) from 2=0.35 to 2=1.4, suggesting that IR galaxies 
dominate Hir at higher redshift. We estimated that Hir includes 
85% of the cosmic star formation at redshifts less than 1, and 
more than 90% at higher redshifts. These r esults are consistent 
with our previous work dCoto et alj|2010bl) . but with much re¬ 
duced statistical errors thanks to the large area coverage of the 
CFHT and AKARI data. 
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